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SOLVING THE KORTEWEG-DE VRIES EQUATION BY ITS 
BILINEAR FORM: WRONSKIAN SOLUTIONS 

WEN-XIU MA AND YUNCHENG YOU 



Abstract. A broad set of sufficient conditions consisting of systems of linear 
partial differential equations is presented which guarantees that the Wron- 
skian determinant solves the Kortcwcg-de Vries equation in the bilinear form. 
A systematical analysis is made for solving the resultant linear systems of 
second-order and third-order partial differential equations, along with solution 
formulas for their representative systems. The key technique is to apply varia- 
tion of parameters in solving the involved non-homogeneous partial differential 
equations. The obtained solution formulas provide us with a comprehensive 
approach to construct the existing solutions and many new solutions including 
rational solutions, solitons, positons, negatons, breathers, complexitons and 
interaction solutions of the Korteweg-de Vries equation. 



1. Introduction 

Among integrable equations is the celebrated Korteweg-de Vries (KdV) equa- 
tion, which serves as a model equation governing weakly nonlinear long waves whose 
phase speed attains a simple maximum for waves of infinite length ,DJ| . It motivates 
us to explore beauty hidden in nonlinear differential (and difference) equations. The 
remarkable and exceptional discovery of the inverse scattering transform |ASe| is 
one of important developments in the field of applied mathematics, which comes 
from the study of the KdV equation. There are various algebraic and geomet- 
ric characteristics that the KdV equation possesses, for example, infinitely many 
symmetries and infinitely many conserved densities MGK , the Lax representation 
|L"]. bi-Hamiltonian structure |Mag| , loop group |S W| . and the Darboux-Backhmd 
transformation MS . More significantly, many physically important solutions to the 
KdV equation can be presented explicitly through a simple, specific form, called the 
Hirota bilinear form. Such exact solutions contain solitons H , S , rational solutions 
[Alga], positons iMatl IR]. negatons |RSK| . breathers [jj and complexitons 

Let us consider the KdV equation in its standard form 
(1.1) u t - 6uu x + u xxx = 0. 

Hirota H introduced the transformation (here and in the rest of the paper, we use 
the notation In / for simplicity, which does not imply that / > 0) 

2(//x* - fl) 



= -2^1n/ = -2(ln/) ra = 



P 
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and transformed the KdV equation into the bilinear form 

(1.2) (D x D t + l>],f ■ f = /,,/ - ft f x + /,,,,/ - i. /•,,,/, + 3/^ = 0, 

which is called the bilinear KdV equation. Here D x and D t are Hirota's bilinear 
operators |BC| . defined by 

f(x + h,t + k)g(x-h,t-k)= V —ADIDlf ■ g)h}k\ 

e - — ' i\v. 

i,j=0 J 

which are also closely related to the vertex operators from elementary particle 
theory. Strictly speaking, we have 

-2(D x D t + D x )f-fi 



u t — 6uu x + u x 



P 



Therefore, if / solves the bilinear KdV equation 1)1.2(1 . then u = —2d x \nf solves 
the original KdV equation (|l.l|l . The bilinear form (|1.2|l looks a little bit more 
complicated than the KdV equation l|l.l|) itself, but its bilinear property, the nearest 
neighbor of the linear property, brings us great convenience in constructing explicit 
exact solutions. Such a bilinear form exists not only for the KdV equation, but 
also for many other physically important soliton equations, such as the Ablowitz- 
Kaup-Newell-Segur equations AKNSj, the Kadomtsev-Petviashvili equation |KPj . 
and the Benney-Roskes equation BR] (or the Davey-Stewartson equation |USj ). 

The Wronskian technique is a powerful tool to construct exact solutions to bi- 
linear differential (and difference) equations. To use this technique, we adopt the 
compact Freeman and Nimmo's notation |FN| : 
(1.3) 



(V -1) = (N- 1;$) =W{<h,<h, ■ 



4°) 

4 0) 



4 (0) 
} N 



4 1 ' 



6 (1) 



'N 



N > 1, 



where 



(1.4) 



j > 1, 1 < i < N. 



Satsuma, Freeman and Nimmo |5J IFN) proved that multi-soliton solutions to the 
KdV equation can be expressed through the above Wronskian determinant, and 
afterwards, Matveev |Mat) generalized the Wronskian determinant which allows 
to present another important class of exact solutions, called positons, to the KdV 
equation. In using the Wronskian method to solve the KdV equation, one usually 
starts from 



(1.5) 



= Xi4>i, <t>i,t = - z 



1 < i < N, 



where Xi are arbitrary real constants, to guarantee that the Wronskian determinant 
(|1.3(l solves the bilinear KdV equation i|1.2fl and thus generates a solution 

(1.6) u = -2fl£lnWfai, fa- ■■ ,M 

to the KdV equation l|l.ljl . The solution determined by a Wronskian determinant 
is called a Wronskian solution. For example, / = W(<fii,<j)2,- • ■ ,4> n ) an d u = 
—2d x h\W(4>\, <f)2, • • • , 4>n) are Wronskian solutions to the bilinear KdV equation 
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1)1.2)1 and the KdV equation , respectively, provided that 4>i , <f>2, • • • ,<f>n satisfy 
1)1. Sirianunpiboon et al. |SHR extended the conditions 1)1.5)1 to the following 

i 

(1-7) -<j>i,xx = Ajj^j, (j)i t t = -^i^xxx, 1 < i < N, 

i=i 

where A^ are arbitrary real constants, in order that the Wronskian determinant 
can generate rational function solutions and their interaction solutions with multi- 
solitons. The first group of the above conditions 1)1.7)1 is a triangular system of 
ordinary differential equations. Thus, solving the system of differential equations 
()1.7)) goes exactly in the same ways for scalar differential equations. 

In this paper, we would like to construct explicit exact solutions to the KdV 
equation 1)1.1)1 by its bilinear form and to illustrate the entire process of construction 
with a broader set of sufficient conditions on the Wronskian solutions: 

N 

(1-8) ~<t>i,xx = X! Aii^j, <f>i,t = ~Hi,xxx + 1 < i < N, 

i=i 

where Ay and £ are arbitrary real constants. Here the first group of conditions is 
a coupled system of ordinary differential equations of second-order, and thus needs 
a more general consideration for solving the involved differential equations. Note 
that the coefficient matrix A = (Ay ) in 1)1.7)1 has only real eigenvalues, but the 
coefficient matrix A = (Ay) in 1)1.8)1 in general has complex eigenvalues. This is 
a significant difference between two sets of conditions, defined by 1)1.7)1 and 1)1.8)1 . 
The substantial extension in ()1.8)) allows us to get a much broader class of explicit 
exact solutions to the KdV equation 1)1.1)1 . In particular, new exact solutions called 
complexitons can be generated from the set of conditions 1)1.811 , but not 1)1.7)1 . 

Now the whole problem of constructing the Wronskian solutions to the KdV 
equation 1)1.1)1 reduces to solve the coupled system of partial differential equations 
l)1.8)l . which turns out to be an interesting mathematical problem itself. We shall 
analyze solution structures on the system 1)1.8)1 in detail, to obtain solution formulas 
for all cases of the system (11.811 . The resultant fundamental solution formulas 
provide a direct and comprehensive approach to construct diverse exact solutions 
to the KdV equation (|1 . II) . such as rational solutions, solitons, positons, negatons, 
breathers, complexitons, and more generally their interaction solutions. 

2. Sufficient conditions on Wronskian solutions 

We begin with stating a broad set of sufficient conditions which make the Wron- 
skian determinant a solution to the bilinear KdV equation 1)1.2)1 . 

Theorem 2.1. Assume that a group of functions <pi = </>j(x, i), 1 < i < N, satisfies 
the following two sets of conditions 

N 

(2.1) -<Pi,xx = ^2Xij{t)<Pj, l<i<N, 

(2.2) fa = -Afa xxx + £(*)&, 1 < i < N, 

simultaneously, where \ij(t) are arbitrary differentiable real functions oft and 

is an arbitrary continuous real function oft. Then f = (N — 1) defined by il.Sp 
solves the bilinear KdV equation 
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The conditions (|2.1|l and Ij2.2(l are a generalization to <|1.7|) presented by Siria- 
nunpiboon ct al. SHR . Actually, if £ = and the coefficient matrix A = (Ay ) 
of (|2.1|l is lower-triangular and independent of time t, then the result of the above 
theorem boils down to the result in SHR . The proof of the theorem needs some 
basic equalities, and the following notation is helpful in our deduction and analysis: 

(N~^~f^ Mi,-" = |$ (0) ,$ (1) ,--- ^W-j- 1 ), ,$(*<)| 

(2.3) = det($ (0) ,$ (1) ,--- ^(^-J- 1 ),^^),... 1 < j < N- 1, 

where ii,--- ,ij are non-negative integers, and the vectors of functions are 
defined by 

(2.4) = (#,#,• • • , <^) T , < j < A 1. 

We also use the assumption for convenience that if i < 0, the column vector $W 
will disappear in the determinant det(- • • , $W ) . . . ). 

Lemma 2.2. Under the conditions ^.lp . the following equalities hold: 

N 

(2.5) AiiW(iV^l) = (jV^3, -/V — 1, A) — (A - 2, A + 1), 



=i 



d 2 (N~^~i) N - - 

(2.6) K — 2 ' = - J2 Ai<(*)(JV - 1) + 2(A — 3, A — 1, A), 



JV 



(2.7) x ti{t)( N - 3, A - 1, A) = (A - 5, A - 3, A - 2, A - 1, A) 

i=l 

+ (A~^3, A, A + 1) - (A~^3, A - 1, A + 2), 



/v 



and 

(2.8) Aii(i)) "(A^l) = (A^5, A - 3, A - 2, A - 1, A) 

i=l 

-(A^4, A - 2, A - 1, A + 1) - (A~^3, A - 1, A + 2) 
+2(A~^3, A, A + 1) + (A~^2, A + 3). 
Proof. Note that we have 

fe=i fc=i »,i=i 

where A = {aij)^xN> an d and Aij denote the determinant resulting 

from \A\ with its fcth row differentiated I times with respect to x, the determinant 
resulting from \A\ with its fcth column differentiated I times with respect to x, and 
the co-factor of a^j, respectively. Choose \A\ as (A — 1) and (A — 3, A — 1, A), 
and use the above equality with I — 2 and the conditions 1)2. ljl . Then we obtain 
the required equalities 1|2. 5JI and 1)2. 7|) immediately. A combination of the equality 
(|2.5(l and the equality 

y gx2 1 = (A - 3, A - 1, A) + (A - 2, A + 1) 
leads to the equality (|2.6I) . 
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We now differentiate (|2.5|) twice with respect to x and utilize 1)2. 6|1 and 1)2.7)1 . 
Then a further computation yields the equality 1)2.8)1 . This completes the proof of 
the lemma. □ 

Proof of Theorem \2.1\ By using the conditions 12. 2|) . we obtain that 

ft = N£f - 4[(JV^4, N - 2, N - 1, N) 

-(iV— 3, N-1,N+1) + (N^~2, N + 2)], 

/x* = N£f s - 4[(jV^5, N - 3, N - 2, N - 1, iV) 

-(iV^3, AT, TV + 1) + (A^2, iV + 3)]. 

Therefore, we can further deduce that 

A := D x (D t + Dl)f ■ f = f xt f ~ f t f x + /,,,,/ - I/,,,/, + :!./;., 

= [-3(jV^5, N — 3, N — 2, N — 1, N) + 6{N~^3, N, N + 1) - 3(7V~^2, TV + 3) 

+3(A^4, - 2, AT - 1, N + 1) + 3(W~^3, N-1,N+ 2)](N~^1) 

-12(]V^3, iV - 1, iV + l)(iV^2, A 7 ") 

+3[(iV^3, N — 1, N) + (iV^2, N + l)] 2 . 

Now using the equalities (|2.5|) and 12.8|l in Lemma \'2. 21 and the Laplace expansion 
of determinants about the last A^ rows, we have 

A = -12(N~^3, N —1,N+ l)(jV^2, AT) 

+12(A^3, N, N + l)(jV^l) + 12(jV^3, N - 1, N)(N~^2, AT + 1) 
A^3 AT -2 AT-1 TV AT + 1 
7V~^3 Af-2 AT-1 N N+l 
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0. 



This shows that / = (JV — 1) solves (|1.2ll . Note that the entries in the last abbrevi- 
ated 2A^ x 2A" determinant above are derivatives of e.g., AT — 1 denotes $( A ' _1 ). 
The proof is finished. □ 
Theorem 12.11 tells us that if a group of functions tfii(x,t), 1 < i < N, satisfies 
the conditions (|2.1|) and 1)2.2)1 . then we can get a solution / = W(<fci, (f>2,- • • , 4>n) 
to the bilinear KdV equation (|1.2fl . The conditions l|2.2|) and (|2.1|) are two linear 
systems of second-order and third-order partial differential equations. We are going 
to solve these linear systems explicitly. It is rather difficult to deal with the case 
where the coefficient matrix A of (|2.1ll has complex eigenvalues. However, when — A 
is a positive definite constant matrix, the situation is easy and the solution formula 
of (|2.1II and (|2.2|) can be expressed as 

(2.9) $ = e-fo £(*') dt' e V^Ax+4A v / ^A* $0j 

where y— A is the square root of —A and $o is an arbitrary initial vector. The 
corresponding Wronskian solutions to the KdV equation , 

(2.10) u = -2d 2 In/ = -2<9 2 ln(A^l;$), 

contain soliton and negaton solutions. 

Before we proceed to solve (12.1)1 and 1)2. 2[l . let us observe the Wronskian deter- 
minants and solutions more carefully. 
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Observation I: From the compatibility conditions (f>i t txx = 4>%xxt, I < i < N, of 
the conditions (|2.1(l and Ij2.2(l . we have the equality 

N 

(2.11) J2 Xi ^ = °' 1 - 1 - N > 

and thus it is easy to see that the Wronskian determinant W(<f>i, 4>2, ■ ■ ■ , 4>n) be- 
comes zero if there is at least one entry Ay satisfying Xij.t 7^ 0. 
Observation II: If we make use of the transformation 

(2.12) j =e-J"o«*)*0 jl l<i<N, 
the differential equations in l|2.2J) can be put into 

<t>i,t = -4<f>i,xxx, 1 < i < N, 

but the differential equations (|2.1|l do not change. Obviously, the resultant Wron- 
skian solutions to the KdV equation are the same: 

(2.13) u=-2d*lnW{<f>i,<f>2,--- ,M = u=-2dllnW{4>!,4>2,-- - ,<M- 

But the Wronskian determinants W((j>i,(j)2, - ■ • ,4>n) and W(0i, 02, ■ ■ ■ , 0/v) are 
different, and thus they gives rise to different solutions to the bilinear KdV equation 

O- 

Observation III: If the coefficient matrix A = (Ay) is similar to another matrix 
M = {nij) under an invertible constant matrix P, let us say A = P~ l MP, then 
<l = P$ solves 

(2.14) $> xx = M$, $ t = -4$^ + 

and the resultant Wronskian solutions to the KdV equation are also the same: 

(2.15) u(A) = -2<9| In |$(°), ,$( Ar - 1 )| 

= -2^1n|P$(°),P$( 1 ),... ,P& N -V\=u(M). 

It follows from Observations II and III that different Wronskian determinants 
may lead to the same solution to the KdV equation. Therefore, based on Observa- 
tion I, in order to construct Wronskian solutions to the KdV equation by its bilinear 
form, we only need to consider the reduced case of l|2.2|l and (|2.1|) under £ = and 
dA/dt — 0, i.e., the following conditions 

iV 

(2.16) ~<t>i,xx = X! ^ijfa' = -40i,axx, 1 < » < N, 
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where Ay are arbitrary real constants. On the other hand, the Jordan form of a 
real matrix has the following two types of blocks: 



(2.17) 



(2.18) 



1 A, 


M 

h Ai 








1 A,- 



ki xki 



h At 



A, 



on 


-Pi ' 




' 1 


" 






, h = 






_ A 









1 



where Ai, and (3i > are all real constants. The first type of blocks has the real 
eigenvalue A^ with algebraic multiplicity fcj, and the second type of blocks has the 
complex eigenvalues \f = cti± [3i\/—l with algebraic multiplicity U. Note that an 
eigenvalue of the coefficient matrix A = (Ay) is also an eigenvalue of the Schrodinger 
operator —-§^2 + u with zero potential. We will present solution formulas for the 
system of differential equations defined by l|2.16|l . according to the situations of 
eigenvalues of the coefficient matrix. Now, based on Observation III, all we need to 
do is to solve a group of subsystems of the sufficient conditions on the Wronskian 
solutions, whose coefficient matrices are of the forms (|2.17|) and i|2.18[l . 
Let us summarize the above analysis as the following theorem. 

Theorem 2.3. For the KdV equation M.l)) . all Wronskian solutions generated from 
the conditions )2.1)) and \2.2]) are among the Wronskian solutions associated with 
the special cases of the conditions $2.16)) whose coefficient matrices A = (Ay) are 
of the Jordan form consisting of two types of Jordan blocks in 12.17]) and 12. Iffy . 

In the next section, we shall consider how to solve the two types of subsystems 
associated with the Jordan blocks in (|2.17|) and l|2.18|) . and present solution for- 
mulas for their representative systems generating Wronskian solutions to the KdV 
equation I|1.1J) . 



3. Solution formulas for the representative systems 

In this section, we would like to consider the construction of solutions to the 
associated system of differential equations defined by (|2.16|l . Based on the form of 
two types of Jordan blocks in (|2.17|) and (|2.18|) . a basic idea to solve the system 
(I2.16|l is to use a recursion process. That is to solve (|2.16|) from <fii to 0jv, one by 
one when the coefficient matrix A = (Ay) has real eigenvalues or pair by pair when 
the coefficient matrix A = (Ay) has complex eigenvalues. Therefore, the entire 
problem is divided into two subproblems - to solve the following two representative 
systems of non- homogeneous differential equations: 



(3.1) 



-<f>xx =X(f> + f, <f>t = 
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and 

!—<h-,xx = Q<01 - (3(1)2 + fi, 
-h,xx = 0<h + a<f> 2 + / 2 , 
= -40i )Xa;x , i = 1, 2, 

where A, a and /3 > are real constants, and / = f{x,t), fx = fx(x,t) and f% = 
fa{x, t) are three given functions satisfying the compatibility condition gt — — ^Qxxx- 
The whole system (|2.16j) whose coefficient matrix is of Jordan form decouples into 
a group of systems determined by (|3.1|) and (|3.2|1 . 



3.1. The case of real eigenvalues. First, let us consider the first representative 
system (13.1(1 . In terms of the eigenvalue A, we will establish solution formulas for 
three situations of the representative system (|3.1(l . 

Zero eigenvalue: The representative system (13.111 in this case is 

(3.3) 4>xx = f, <f>t = ~4(j) x xx, 

where / = f(x,t) is a given function satisfying the compatibility condition f t — 
—'If xxx- Immediately from the first equation in (|3.3|) . we can get 

f-X r-x' 

(3.4) <j> = / f(x",t)dx"dx' + c(t)x + d(t). 

Jo Jo 

Note that ft = — 4/ xxx , and thus the second equation of (|3.3|) equivalently requires 

Ct = — 4:fxx\x = 0, <lt — —4:fx\x=0- 

Then, we obtain 

(3.5) c{t) = -if f X x(0,t')dt' + c , d(*) = -4 / f x (0,t')dt' + d , 

Jo Jo 

where cq and do are arbitrary real constants. Summing up, we have the following 
theorem. 

Theorem 3.1. Let f = f(x,t) be a given function satisfying the compatibility 
condition f t = —Af xxx . Then the system of differential equations has the 

general solution given by {3.$ and f,y.,5J) . 

Negative eigenvalue: In this case, the representative system (|3.1f) reads as 

(3.6) <p xx = a<j) + f, <f> t = -4:(j> xxx , a > 0, 

where a is a constant and / is a given function satisfying the compatibility condition 
ft = — 4/ xxx . From two basic solutions e^ x and e~^* x of 4> xx = a<p, we get a 
general solution of the first differential equation of (|3.6|) : 



(3.7) = \J— T f(x>,t)e-^'dx' + c{t) 

12^/aJo 

-i= [ X f(x',t)e^ x 'dx' + d(t) 
V" Jo 
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by variation of parameters. And then noting that ft = —4f X xx, by a direct but 
lengthy computation, we find that the second differential equation of (|3.6|l equiva- 
lently requires 

c t = -4oty/ac - 2(-^= f xx + f x +v / "/)|x=o, 
V« 

d t = 4a\fad+2{ l — f xx + f x - Va/)|s=o< 

These are two linear systems for c and d, respectively. Hence, we can immediately 
obtain the solution formulas for c and d: 

(3.8) c(t)=e- 4a ^ t \c -2 f e 4a ^*'(-^ f xx + f x + ^f) {Q,t'))dt' 



o 



(3.9) d(t)=e 4a ^ t d + 2 fV 4 ^''!- -L fx X + fx-V^f)M)dt' 

'a 



where cq and do are arbitrary real constants. We conclude this result as follows. 

Theorem 3.2. Let f = f{x,t) be a given function satisfying the compatibility 
condition f t = —Af xxx . Then the system of differential equations \3. 6)) has the 
general solution given by \S. 7[ ) ; and ^3.y\) . 

Positive eigenvalue: The representative system (|3.1ll in this case is 

(3.10) <j> xx =-ouf> + f, a> 0, 

(3.11) 4>t = -4:<t>xXX, 

where a is a constant and / = f(x,t) is a given function satisfying the compati- 
bility condition f t = —4:f xxx . To solve the differential equation l|3.10[) . we can start 
from two basic solutions sin y/ax and cos yfax of the corresponding homogeneous 
equation of (|3.1Q(I . and then use the method of variation of parameters. A general 
solution is listed in the following lemma for reference. 

Lemma 3.3. For any given function g = g{x) and any positive constant a, the 
second-order ordinary differential equation tp xx = —aip + g has the general solution 

f x sin y/a(x — x 1 ) . .. , . ,— 

ip = I -= g(x)dx + c cos v era + d sm v ai, 

Jo V a 

where c and d are arbitrary real constants. 

Now to solve the partial differential equation (|3.11|) . we start with the solution 
of H3.10f) given by Lemma 13.31 and use the method of variation of parameters again. 
Therefore, the solution of l|3.11|l is assumed to be 

(3.12) = / v P ' f(x',t)dx' + c(t) cos y/ax + d(t) sin \fax. 

Jo V" 
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Directly we can compute that 

4>xxx = fx- f(x',t)(acosy/a(x- x'))dx' 



+c(t)ay/a sin \fa~x — d{t)ayfa cos yfax, 



SinVra a £ \ /- ami r i— • i— \ 

= ifxx U=o 7= 1- 4/ x x =o cos \Jax - 4:f\x = oW asm v ax ) 

\la 



-4/a,+ / if (x',t)(a cos */a(x - x'))dx' 

+Ct cos \fax + d t sin y/ax, 

where ft = —ifxxx has been used for the computation of t . Therefore, the differ- 
ential equation l|3.11fl equivalently requires 
4 

(3.13) —^fxx\x=o ~ 4^/1^=0 +d t = -Aay/ac, 
V« 

(3.14) 4f gl \ x= o + ct = 4aVad. 
Then it follows that 

c tt + 16a 3 c + h = 0, 
where the function h = h(t) is defined by 

(3.15) h := l6afxx\x=o - 16a 2 /U=o + &f xt \ x=0 . 

Again by Lemma 13.31 we obtain 

. . f sin4a v / a(i — t') . . _ . _ 

cli = — / = h{t)dt + Co cos 4a Vat + do sin 4a\/at, 

J 4a V a 

where cq and do are arbitrary real constants. Noting that h is given by Ij3.15|l . we 
have 

(3.16) c(t) = f x (0,0) Sm4a ^ t -4 f f x (0,t')co S 4aV^(t-t')dt' 

ay/a J 

r* i 

-4 / (— / xx -^/)(0,t')sm4a v / 5(t-t')^ 
■/o V" 

+Co cos Aayfat + do sin Aay/at, 
and then by 13.14|l . we have 

(3.17) d(t) = f x (0, 0) C0S 4 "/" t + 4 f / a (0, sin 4a^(t - t') dt' 

a\/a Jo 

rt 



r i 

-4 / (-=/„- V^/)(0,t')cos4a v / ^(t-t')*' 
Jo V" 

-Cq sin Aay/at + do cos 4av / a£, 



where Co and <io are independent of x and i. 

Finally summing up, we have the following theorem on the general solution of 
the system of differential equations (|3.1()|l and H3.11JI . 

Theorem 3.4. Let f — f(x,t) be a given function satisfying the compatibility 
condition ft — —Afxxxi and a be a positive constant. Then the system of differential 
equations ^3. lip and $3.1U\) has the general solution given by \3.1'ify with c(t) and 
d(t) being shown by \3.16\) and {&. 1 7| ), respectively. 
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3.2. The case of complex eigenvalues. Let us now consider the second repre- 
sentative system (13.2(1 . To present its solution formula, we first solve the following 
coupled system of second-order, non-homogeneous ordinary differential equations 



(3.18) 



h,xx = ot<t>i ~ + fl, 

t>2,xx = P4>1 + a<j>2 + /2, 



where a and (3 > are real constants, and f\ and f 2 are two given functions of x. 
Since (3 > 0, its coefficient matrix 



.4 



a -P 
P a 

has two complex eigenvalues X^ = a ± (3i. 

Theorem 3.5. The coupled non-homogenous system of second- order ordinary dif- 
ferential equations has the general solution 



(3.19) fa 

(3.20) ^ 

(3.21) 4> h 2 



(Di cos 5x + D 2 sin 5x) e Ax + (D 3 cos Sx + D4 sin Sx) e Ax , 
(-D 2 cosSx + Di sm5x)e Ax + (— -D4 cosSx + D 3 sinSx)e~ Ax 
1 



[fi{y)cosS(x ~y)~ f 2 (y) sinS(x - y)] sinh A(x - y) dy, 

1 f x 



(3.22) # 

(3.23) # 
where two constants A and 5 are given by 



(3.24) 



\J a 2 + [3 2 — a 



^a 2 + I3 2 + a 



and Di, 1 < i < 4, are arbitrary real constants. 



Proof. First let us consider the homogenous case of l|3~T3jl . i.e., the case of ((3.18(1 
with /1 = /a = 0. Therefore, we have 



(3.25) 



1 . <* . 

h = S<Pl,xx + -£<Pl- 

p p 



Then replacing 4> 2 with this expression in the second equation of the homogenous 
case leads to the following fourth-order differential equation 



0. 



(3.26) 4>i,xxxx + 2a4>i,xx + (< 

Its characteristic polynomial 

p(\) = A 4 + 2aX 2 + (a 2 + (3 2 ) = 

has four complex roots 

Xf = A± SV-L, Xf = -A ± SV^l, 

where A and S are defined by (|3~2l|) . Therefore, we have the general solution <f>i 
and (^2, as given in the theorem, for the homogeneous case of 1(3.18(1 . 
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To construct a special solution, 0f and <f> 2 , for the non- homogenous case of l|3.18jl . 
we adopt the method of variation of parameters for Di, 1 < i < 4, and thus assume 
that 

{<j)\ = (DxCOsSx + D 2 sin5x)e Ax + (D 3 cosSx + D 4 sin Sx)e~ Ax , 
A -A 
<j>2 — (— D2 cos 5x + D\ sinfa)e + (— D 4 cos fa + D 3 sinfa)e , 

where Di, 1 < i < 4, are viewed as functions of a;. Then we can make the following 
computation 

(3.28) (f> s hx = [(AD 1 +5D 2 )cos5x+(AD 2 -8D 1 )sm5x]e Ax 

+ (5D 4 - AL> 3 ) cos 5x - {SD 3 + AD 4 ) sin fa] e" Ax , 

(3.29) <f>i tXX = [A(AD 1 +6D 2 ) + 6(AD 2 -SD 1 )](cosSx)e Ax 

+ [A(AD 2 -6D 1 )- S{AD 4 + 5D 2 )} (sin 6x)e Ax 



[S(6D 3 + AD 4 ) - A(SD 4 - AD 3 )] (cos 5x) 



-Ax 



+ [S(5D 4 - AD 3 ) + A{SD 3 + AD 4 )] (smSx)e- Ax 
= {aD 1 + (3D 2 )(cos5x)e Ax + (~(3D 1 + aD 2 ){smSx)e Ax 

+ (aD 3 - PD 4 )(cosSx)e~ Ax + {f3D 3 + aD 4 ){sin5x)e- Ax 
:= (Di cosSx + D 2 sindx)e Ax + (D 3 cos fa + D 4 sin 6x)e~ Ax , 

(3.30) <$> s %x = [{-AD 2 + 6D 1 )cos5x + (AD 1 +5D 2 )siii6x]e Ax 

+ (<5 D 3 + AD 4 ) cos fa + {SD 4 - AD 3 ) sinfa]e _Ax , 

(3.31) <t> s 2tXX = {-aD 2 +f3D 1 )(cos5x)e Ax + {(3D 2 +aD 1 ){sm5x)e Ax 

-{aD 4 + f3D 3 )(cos6x)e- Ax - {J3D 4 - aD 3 ) (sin Sx)e~ Ax . 
In the above computation, the method of variation of parameters required that 
{D\ lX cos fa + D 2 . x sin Sx)e Ax + (D 3 . x cos Sx + D 4 , x sin Sx)e~ Ax = 0, 
[(AD lyX + 5D 2 ^ X ) cos fa + (A£>2,x - SDx^) sinfa]e Aa: 
+ (5D 4 . X - AD 3>X ) cos fa - (SD 3:X + AD 4 , X ) sin5x]e~ Ax = - fx, 
{—D2, x cos fa + D\ lX sin 8x)e Ax — (D 4 . x cos fa — D 3 . x sin fa)e~ Ax = 0, 
[(-AD 2>X + SDx >x )cos6x + (AD lyX + 6D 2 . x )smSx]e Ax 
+ (5D 3 , X + AD 4 , X ) cos fa + (SD 4iX - AD 3iX ) sin fa] e - Ax = -fa- 
Solving this linear system for D^ x leads to the expressions: 

Di, x = - — (frcosSx + f 2 sin5x)e~ Ax , 
D 2 . x = sinSx - f 2 cos5x)e~ Ax , 

D 3:X = -^-(/i cos fa + f 2 sm5x)e~ Ax = -Di. x , 



2A 

D 4 . x = -^-(/lsinfa - f 2 cosfa)c~ A ' T = —D 



2A 

Then integrating them to obtain the expressions for Di and inserting the resultant 
expressions of Di into Ij3.27|l yield the special solution, 0f and <p 2 , given in the 
theorem. The proof is finished. □ 
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Let us then present the solution formula for the second representative system 
(|3.2|1 . which can be used to generate new type Wronskian solutions to the KdV 
equation. 

Theorem 3.6. Assume that f\ = f\(x,t) and / 2 = f2(x,t) are two given functions 
satisfying the compatibility conditions f± t t — —^f\. xxx and / 2jt = —^fi,xxx- Then 
the general solution to the system of differential equations is given by $3.19\) - 

\S. 2<il\) with the coefficients Di = Di(t), 1 < i < 4, being determined by 



(3.32) 





= e &t 









cos [3t — sin [3t 
sin (3s cos (3t 



D 20 



cos [3(t — s) — sin (3(t — s) 
sin/3(t — s) cos (3(t — s) 



(3.33) 



£3 " 


= e~ at 


cos fit 


— sin (3t 




' D 3Q ' 


D 4 _ 




sin (3s 


cos (3t 







-a(i-s) 



cos (3(t — s) — sin (3{t — s) 
sin/3(t — s) cos(3{t — s) 



P2(s) 



Ps(s) 



ds. 



ds. 



where D^, 1 < i < 4, are arbitrary real constants, and two constants a and (3 are 
given by 

(3.34) a = -4A(Va 2 + (3 2 + 2a), (3 = A6(y/a 2 + (3 2 ~ 2a), 

and four functions, pi, 1 < i < 4, are defined by 



2 2a 2S 

-^fl,xx\x=0 + 2fi tX \ x= Q + -^-/l|a;=0 — ~^f2,x\ x =0 ~ 4(5/ 2 U=0; 



-^fl,x\x=0 ~ 4<5/l|o;=0 — -^f2, XX \ X =0 — 2/ 2 ,xU=0 ~^f2\ X =0: 



fl,xx\x=0 + 2/l ; x|x=0 



2a . 2(5 

-^"/l|x=0 + "^"j2,a:U=0 — 4 0,/2|x=0, 



Pl(t) 
P2(t) 

Pa(*) 

P4(*) 

wif/i A and 5 fcinj defined by 

Proof. Note that now /1 = /i(a;, i) and / 2 = /2(^, are functions of two variables 
x and t. But in the solution formula (|3.22|) and l|3.23[l . t is the dummy variable. 
Let us first compute that 



25 2 2a 

~^fl,x\x=0 — 45/l|x=0 + -^f2,xx\x=0 — 2f2,x\x=0 + ~^f2\x=0, 



b\ t = (Di^t cos Sx + D2,t sin Sx)e + (Ds t t cos Sx + D^t sin Sx)e 

!■!,< = - [ [9i(x-y)fi.t(y^) + 92(x-y)f2,t(y^)] d y 

[gi(% ~ y)fl,yyy{y,t) + g 2 (x - y)f2,yyy{y 7 t)] dy 



-Ax 



-4fll(a:)(/i,xx|x=o) 

~452(a:)(/2,a;a:U=o) 



7i,x(ar)(/i,xU=o) 
?2,x(a;)(/2,x|x=o) 



?i,xx(a;)(/i|x=o) 
?2,xx(a;)(/2|x=o) - Wi, xxx , 
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<7i = — (cos <5x)(sinh Ax), g 2 = — ^(sin^^) (sinh Ax). 



where 



Second, based on (|3.28|> and (|3.29|) . we have 

§\, xxx = [(ADi + SD 2 ) cosdx + {AD 2 - 5Di) sinSx]e Ax 

+ [(5Di - AD 3 )cosSx - (SD 3 + AD 4 ) sinfe]e- Ax , 
and a direct computation gives rise to 

9iA x ) = l(™s6x)(e Ax + e~ Ax ) - ^ (sin Sx) (e Ax - e~ Ax ), 

gi,xx{x) — {— SwaSx + cosSx)e Ax — (Ssin5x + -—^ cosSx)e~ Ax , 



92,x( X ) 



1 



(sinfeXe^ + e"^) 



-Ax 



2A 



(cos Sx) (e 



Ai e _Aa 



92,xx(x) — —(ScosSx + sin<Sx)e Ax — (S cosSx — sin5x)e Ax . 
Now a comparison of coefficients of four functionally independent terms 
(cos5x)e Ax , (sinSx)e Ax , (cos5x)e~~ Ax , (sinSx)e~ Ax 



in the equation 



n,t — 



i.e. 



•) n,t -r vi,t - —*4>i, XX x ~ ^Ixxxi yields that 

-4(ADi+<y^ 2 )+Pi(*) 

4(3A<5 2 - A 3 )Z?i - 4(3A 2 <S - <5 3 )L> 2 + pi(t) 

a£>i -0D a +Pi(t), 

-4(AL) 2 -^ 1 )+p 2 (t) 

4(3A 2 <S - <S 3 )L>i + 4(3A5 2 - A 3 )L> 2 + p 2 (t) 

pD 1 +&D 2 +p 2 {t), 

-4(5D A -AD 3 )+p 3 (t) 

-4(3A6 2 - A 3 )L> 3 - 4(3A 2 (5 - S 3 )D 4 + p 3 (t) 

-&D 3 -pD 4 +p 3 (t), 

-4(-SD 3 -AD A )+p 4 {t) 

4(3A 2 6 - S 3 )D 3 ~ 4(3A<5 2 - A 3 )L> 4 + p 4 (t) 

)9D3-aD4+p4(i). 



(3.35) 



(3.36) 



(3.37) 



(3.38) 



D 



i,t — 



Do 



D 



3,t — 



£»4 



These are two linear systems of first-order ordinary differential equations. One is 
for Di and D 2 and the other is for D 3 and D4. Solving these two coupled linear 
systems for Di, 1 < i < 4, we obtain the general solution given by (|3.32|) and (|3.33() 
in the theorem. For example, we rewrite (|3.35|) and 1)3. 36[) as 







' D 1 ' 




Pi 




~ a -0 ' 


= A 




+ 


, A = 




. D U . 








. P 2 . 




. & 



Then it follows that 



" Di ' 


= e At 


' D w ' 


+ e At [ e~ As 


Pi( s ) 








Jo 





ds . 
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Since A has two complex eigenvalues 

A ± =a±/3 v /Z T, 

the evolution operator e At reads as |HS1 IAP| 

. . T cos Bt — sin Bt 
e At = e at _ 

sin Bt cos Bt 

The derivation of the formulas for D 3 and D4 is completely similar. Note that the 
other equation 2 ,t = —^4>2,xxx is automatically satisfied, once Z^, 1 < i < 4, are 
determined as above. The proof is finished. □ 

4. WRONSKIAN SOLUTIONS 

In principle, we can construct general solutions associated with two types of Jor- 
dan blocks of the coefficient matrix through the solution formulas established in the 
previous section. In what follows, we present a few special but interesting Wron- 
skian solutions, rather than general ones which are too complicated and lengthy to 
present. 

4.1. Rational solutions. Let us start with 

(4.1) ip( v ) = e"*- 4 " 3 * - e -"* +4 " 3t = 2sinh (r}x - 4ry 3 t), 

where 77 is a real constant. It is easy to sec that this function t/jiv) solves 



(4.2) 

Upon expanding 
(4.3) 



Ipxx =Ti}, 1p t = ~^xxx- 



V(r/)=]T<M 2l+1 , 

i=0 



from l|4.2ll we obtain 

(4.4) cf> 0tXX = 0, 4> i+ i,xx = <pi, 4>i,t = ~^4>i,xxx, i > 0. 

The expressions of tfii can be calculated through the series expansion of the hyper- 
bolic function sinhx: 

[(2i+l)/3] 



(4.5) 



E 

3=0 



(-Ay 



& - 3j + iy.j\ 



c 2i - 3j+1 t j , 1 > 0. 



Then for each k > 1, the associated Wronskian solution to the KdV equation (|l.l(l 
is given by 

(4.6) u= -2dl In W(<f> , fa,--. 

Obviously, this solution is rational and it corresponds to the following Jordan 
block: 



(4.7) 





1 







1 



fexfe 
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We call it a rational solution of order k— 1. This type of solutions was also discussed 
in jASal ISHRI IAM| and obtained by Freeman and Nimmo |FNj by taking the limit 
r/i — > of the multi-soliton solution. Two rational solutions of lower-order are 



-2d 2 x \nW(4> ) = - 
-2^1nW(^,,0i) 



2d*lnW(x) = ^, 

= -2dl\nW(x, -x 3 
6 



At) 



6x 4 - lUxt 



(x 3 + 12t) 2 ' 

4.2. Solitons, positons and negatons. For a nonzero real eigenvalue Ai, we start 
from the eigenfunction 4>i(\i) determined by 

(4.8) -{(f)i{\i))xx = Xi4>i(\i), (&(A<))t = -4(<^ i (A j ))xxx. 
General solutions to this system in two cases of Ai > and Ai < read as 

(4.9) (f>i(Xi) = C u s\n{qix + Arjft) + C 2i cos^a; + 4r]ft), r]i = \/a7, 

(4.10) <t>i{\i) = C lt sinhfox - A-qft) + C 2l cosh^x - 4r)ft), i], = \/^\ , 

respectively, where Cu and C 2 i are arbitrary real constants. By an inspection, we 
find that 



Ai 

1 A,- 







1 Ai 



kiXki 



UK) 

±d Xi 4>i{\i) 



and 



(^>(Ai))t 



4(-^^i(Ai)) xxx , 0< j<h-l, 



J 



where dx t denotes the derivative with respect to Ai and fci is an arbitrary non- 
negative integer. Therefore, through this set of eigenfunctions, we obtain a Wron- 
skian solution to the KdV equation 1)1. lfl : 



(4.11) 



1 



1 



^"Vi(Ai)), 



u = -2d 2 x lnW(MK), ffxMK), ' ■ ■ , (fc ,_ 1)r Ai 

which corresponds to the first type of Jordan blocks with a nonzero real eigenvalue. 

When Ai > 0, we get positon solutions |Mat| . and when Ai < 0, we get negaton 
solutions RSKj . A more general positon or negaton can be obtained by combining 
n sets of eigenfunctions associated with different A; > or different Ai < 0: 



(4.12) 



1 



u = -2^1n^ 1 (A 1 ), ^. ^rViCAi); 



; </>n(A„), 



1 



-d^MK)). 

This solution is called an n-position or n-negaton of order (fci — 1, k 2 — 1, • • • ,k n — 
1), to reflect the number of different eigenvalues and the orders of derivatives of 
eigenfunctions. If ki = 1, 1 < i < n, we simply say that it is an n-position or 
rt-negaton. 

An n-soliton solution is a special n-negaton: 
(4.13) u= -2dl\n{<j>i,<l> 2 ,--- ,(f> n ) 
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with 4>i being given by 

(4.14a) <f>i — cosh(r]iX — Arftt + Ji), i odd, 

(4.14b) (f>i — sinh(r]iX — 4?]ft + 7,), i even, 

where < 771 < 772 < • • • < T) n and ji, 1 < i < n, are arbitrary real constants. There 
are other representations of multi-soliton solutions to the KdV equation (see, for 
example, [Hi IBS| ). Two kinds of special positons of order k are 



(4.15) u = -2d*lnW{<f>,d x <t>,--- .^a~V), <P = cos^x + AtjH + 7(77)), 

(4.16) u = -2dl In W(<j>, 9 A 0, ■ ■ • , d£~ V), ^ = s in (^ + 4? / 3i + 7^)) = 

where A > 0, 77 — \f\ and 7 is an arbitrary function of 77. But these two kinds of 
positons are equivalent to each other, due to the existence of the arbitrary function 
7. Similarly, two kinds of special negatons of order k are 

(4.17) u= -2dl\nW{<f>,d x <j>,--- <t> = cosh^a; - Arft + 7(77)), 

(4.18) u= -2^]hW(^,9a</»,--- >^"V), = smh(77x- 477-^ + 7(77)), 

where A < 0, 77 — \J — A and 7 is an arbitrary function of 77. These solutions exhibited 
above were also discussed in a little bit different way in |APPI IMatl QlSK . The 
following are solitons, positons and negatons of lower-order: 

-2?7 2 

u = — 2<9 2 ln(cosh(7ia; - Arft + 7)) — 



cosh {qx — Arj 3 t + 7) 



2?7 2 

-2<9 2 ln(cos(7ia; + 47y 3 7j + 7)) = 



= —2d* ln(sinh(77x — Arft + 7)) = — 



cos 2 (772; + Arft + 7) ' 
2t7 2 



sinh 2 (772; — 4?7 3 t + 7) ' 
u = -2&1 In V7(cosh(77 1 x — Arjft + 71), sinh^a; — Ar]\t + 72)) 

_ 4(?7 2 -T7|)[(r72 - ill) +Vi cosh 20 2 + rj% cosh20i] 
[(772 - 771) cosh(6»i + 02 ) + (?72 + 771) cosh(0i - 2 )] 2 ' 
U = -2<9 2 In W(cos(r]x + Arft + 7), d\ cos(r)x + Aift + 7)) 

167? 2 [2cos 2 (?7a; + 4?7 3 t + 7) + {r\X + l2r)H) sin 2(770; + 4?7 3 t + 7)] 
~~ [2(772; + 1277 3 t) + sin 2(770; + 4?7 3 t-|-7)] 2 

u = -2d\ In W(cosh(r]x — A^t + 7), d\ cosh(772; — 4rj 3 t + 7)) 

16?? 2 [2 cosh 2 (772; - 4t7 3 < + 7) - (772; - 12rj 3 t) sinh 2(772; - Ai] 3 t + 7)] 
[2(772; - 1277 3 /j) + sinh 2(772; - 477 3 t + 7)] 2 ' 

where 77, 7, 77^ and 74 are arbitrary real constants, and 

0i = mx - 4?7 3 i + 7i , i = 1,2. 

Some parts of the graphs of these solutions with 7 = 7^ = are displayed in the 
following Figures 1 and 2. 
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4.3. Complexitons. For the second type of Jordan blocks of the coefficient ma- 
trix, we start from a pair of eigenfunctions <J>j(ai,ft) = (<fin(ai, f3i), fyaiui, ft)) T 
determined by 



(4.19) 

and 
(4.20) 



(f>n{ai, (3i) 





OH -fa 






Pi Oil 



i(oti,0i))t = -4((f>ij(ai,(3i)) xxx , j = 1,2. 



The coefficient matrix Ai has two complex eigenfunctions: ajifty^-T. By Theorem 
13.61 a general solution to the homogenous system of Ij4.19|) and (|4.2()|l is 

1, 



(4.21) 



(4.22) 



[C0s((5j(x - Pit) + Kli)] 



1, 



+|[cos(« i (x - Pit) + K 2i )] e-^+^t)-^ 
i[sin(^(x - Pit) + KM)] e^+^+T" 
-i[sin(^(x - ftt) + km)] e - A *( a: + a **)-r« , 



where Kjj and 7^ are arbitrary real constants and 
(4.23) 



(4.24) a, - + /?? + 8a t , ft = 4^ a? + ft? - 8a 4 

Similarly by an inspection, we can see that 



4\ 



and 



(2;-l)! U a s 1 



1 



h Ai 







h Ai 



liXli 



1 d'i- 1 * 



1 



J 



J 



where d ai denotes the derivative with respect to ai. Taking the derivative with 
respect to ft, we can have 



V.dp^i 



1 a'i-lfr 



-4, 

S 2 ^ 







S 2 At 



1 ati-lfr 



and 



(— ^*i) t = -4(-^$i)~*. < i < h ~ 1. 



il A 
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where 



S 2 = 



-1 

1 



Therefore, we obtain two Wronskian solutions to the KdV equation ll.lfl : 



(4.25) 

and 
(4.26) 



1 



T 



1 



'(k-l)\ a < 
1 



(k-iy. ft 

which correspond to the same Jordan block of the second type but are different. 
Such a general Wronskian solution is 



(4.27) u = -2d*\nW(<f>i 



6k- 1 & • ■ 

'(h-iy. Cl 



where can be either of d ai and dp i . This solution is called an rt-complexiton 
solution of order (1% — 1, 1% — 1, • • • , l n — 1), to reflect the orders of derivatives of 
eigenf unctions with respect to eigenvalues. If ij = 1, 1 < i < n, we simply say that 
it is an n-complexiton. 

In particular, one-complexiton reads as [M] 

(4.28) u = -2dl\nW(4> u ,^ 2 ) 

-\p\ [l + cos(2<5i (x - (3it) + 2k x ) cosh(2Ai(a; + a x t) + 2 7 i)] 

[Ai sin(2(5i(a; - fat) + 2m) + 5i sinh(2Ai(a; + out) + 2 7l )] 2 

4ai/3i sin(25i(x - Pit) + 2ki) sinh(2Ai(a; + ait) + 2ji) 

[Ax sin(2<5i(.T - pit) + 2Ki) + S x sinh(2Ai(x + ait) + 2 7l )] 2 ' 

where a±, Pi > 0, n x — Ku — K12, and 71 = 711 = 712 are arbitrary real constants, 
and Ai, Si, ai, and pi are given by 



A 1 = 



a 1 



■s/aj+pj - ai 



61 = 



sja\ + P'{ + a x 



>a\ + Pl+ Sax, Pi = iJaj + P\ - 8«i 



The subcase of ai = leads to the following solution 

8P1 + 8P1 cos(V2^l x - 4Piy/2~Pi~t + 2k x ) cosh(V2ft"x + 4/?i ^/2p\t + 2 7l ) 

u = ■ ■ ■ ■ 5 — . 

[sm(y/2(3i~x - APi^JWit + 2k x ) + sinh( v / 2^Ta; + 4/3i V / 2^Ti + 2 7i )] 

This solution is associated with purely imaginary eigenvalues of the Schrodinger 
spectral problem with zero potential. In addition, if we fix that 



Kl 



4' 



7i 



5 *(-), Ai 
I a 



-a, Si 



where a and b are arbitrary real constants, then our one-complexiton solution 
boils down to the breather-like or spike- like solution presented in [[]: 

8{(a 2 — b 2 )(b/a) cos^sinh(?7 + p) + 2b 2 [l + sin^cosh(7y +£>)]} 
[cos v — (b/a) sinh(?7 + p)] 2 
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where v = -2bx + 8(3a 2 b - b 3 )t 7 rj = -2ax + 8(a 3 - 3ab 2 )t and p = ln(6/a). Figure 
3 depicts some singularities of the single complexiton solution with k = 7 = 0. 

4.4. Interaction solutions. We are now presenting examples of Wronskian inter- 
action solutions among different kinds of Wronskian solutions to the KdV equation. 
Let us assume that there are two sets of eigenfunctions 

(4.29) <MA),0 2 (A),--- ,&(A); M»),Mv),--- ,Mv) 

associated with two different eigenvalues A and fx, respectively. A Wronskian solu- 
tion 

(4.30) u = -2^W(0i(A),..-,0 fc (A);Vi(M).-" >z(m)) 

is said to be a Wronskian interaction solution between two solutions determined 
by the two sets of eigenfunctions in Ij4.29|l . Of course, we can have more general 
Wronskian interaction solutions among three or more kinds of solutions such as ra- 
tional solutions, positons, solitons, negatons, breathers and complexitons. Roughly 
speaking, it increases the complexities of rational solutions, positons, negatons and 
complexitons, respectively, to add zero, positive, negative or complex eigenvalues 
to the spectrum of the coefficient matrix. 

In what follows, we would like to show a few special Wronskian interaction 
solutions. Let us first choose different sets of eigenfunctions: 

^rational = CX + d, C, d = COnstS.; 

^soiiton = cosh(ryx - Arft + 7), 77, 7 = consts.; 
P ositon = cos^x + Arft + 7), 7], 7 = consts.; 

\/2 v/2 

^complexiton, 1 = COs(— (x - At) + «) COSh(— (x + At) + 7), K, J = COnstS., 

V2 V2 
^compicxiton.2 = sin(—(x - At) + k) sinh( — (x + At) + 7), K,7 = consts. 

Three Wronskian interaction determinants between any two of a rational solution, 
a single soliton and a single positon are 

Vibrational, Soliton) = v(cx + d) sinh^x - Ar/ 3 t + 7) - c cosh(fya; - Arft + 7), 

Vibrational, ^positon) = -1){cX + d) Sm(r]X + Arft + j) - CCOs(l]X + Arft + 7), 

W"(0 so iiton, ^positon) = -r) cosh^x - Arft + 7) sin(?7x + Arft + 7) 
— rj cos(ryx + Arft + 7) sinh(??x — Arft + 7). 
Further the corresponding Wronskian interaction solutions read as 

9*2 1 WU A \ ^ 2 [V 2 (cX + d) 2 +C 2 CO S h 2 ^} 

U rs = -2djlnVF(0rational,0 S oliton) ~ 



[t)(cx + d) sinh£_ — ccosh£ 

-2&1 In W (^rational, 



_ 2r] 2 [rj 2 {cx + d) 2 - c 2 cos 2 £ + ] 

positon ) "f 7 : "77 : Z j Z y7 ■ 

[j]{cx + d) sm^-|- + ccos^+J^ 



9^2, WU A \ 47 ? 2 (C0S11 2 e_ + COS 2 £+) 

-2d^ln W(0 so liton,<Ppositoii) ~ 



(cosh£_ sin£+ + cos£+ sinh£_ 
where 

£± — V x i Arft + 7. 
Some singularities of these solutions are shown in Figures 4 and 5. 
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The following are two Wronskian interaction determinants and solutions involv- 
ing three eigenfunctions. The first Wronskian determinant is 

WX'/'rationaU ^soliton : ^positon) 

= —i] 2 {(cx + d)i](cos^ + sinh£_ — sin£ + cosh£_) — 2ccos£ + cosh£_] 
so that its corresponding Wronskian solution reads as 

— 2d 2 In Vibrational; "^soliton, ^positon) 

4(cx + d)rj 3 (sm£ + cosh£_ + cos£ + sinh£_) 
(cx + d)r)(cos sinh£_ — sin£ + cosh£_) — 2ccos£+ cosh£_ 

2rj 2 [crj cos sinh£_ + 2(cx + d)rjsin£+ sinh£_ — cr; sin£+ cosh£_] 2 



+ 



[(cx + d)rj(cos £+ sinh £_ — sin £ + cosh ) — 2c cos £ + cosh ] 



2 



where 



£± = tjx ± Arft + 7. 

The second Wronskian determinant and its corresponding Wronskian solution are 

W(x, </> C ompl exiton,U ^complexiton,2) 

\[2 \[2 2 2 

= — — xsin2£ + -^-xsinh2C — cosh ( + sin £, 

and 

u rc 2^ In VT(x, (^» com p} cx iton,l ' ^complcxiton^) 

= [cosh 4 C - cos 4 £ - Ax 2 cos 2 £ cosh 2 C + [2x 2 - f ) cosh 2 ( 

+ (2x 2 + 1 ) cos 2 £ + 4 V2 x cosh 2 C sin 2£ + 4\/2 x cos 2 £ sinh 2( 

1 / /2 

- V2 x sinh 2( - V2 x sin 2£ - - sin 2£ sinh 2Q] j (- -^x sin 2£ 



/2 

+— xsinh2C - cosh 2 Q + sin 2 £) 2 , 



where 



£ = i V2x - 2V2t + k, (=^V2x + 2V2t + 
Figure 6 depicts some singularities of these two solutions. 



5. Concluding remarks 

A broad set of sufficient conditions consisting of coupled systems of linear par- 
tial differential equations has been presented, which guarantees that the Wronskian 
determinant solves the Korteweg-de Vries equation in the bilinear form. A system- 
atical analysis has been made for solving the resultant coupled systems of linear 
partial differential equations, and solution formulas for their representative systems 
have been explicitly presented. The key technique is to apply the variation of pa- 
rameters in solving the involved non-homogeneous partial differential equations of 
second-order and third-order. The whole analysis also offers an approach to solve 
the coupled system of partial differential equations: 

N 

~4>i,xx = ^ *y0J> = ~^4>i,xxx, 1 <i < N, 

1=1 
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where Ay are arbitrary real constants. 

Moreover, for each type of Jordan blocks of the coefficient matrix A = (Ay), 
special sets of eigenfunctions have been constructed and used to generate rational 
solutions, solitons, positons, negatons, breathers, complexitons and their interac- 
tion solutions to the Korteweg-de Vries equation. Of course, the obtained solution 
formulas of the representative systems allow us to construct more general Wron- 
skian solutions than rational solutions, positons, negatons, complexitons and their 
interaction solutions presented in this paper. The resultant solution structures also 
show us the rich diversity that the KdV equation carries. We believe that any 
(2 + l)-dimensional counterpart of the KdV equation will have much more gen- 
eral solution structures. This is because higher dimensional equations have bigger 
spaces of initial data to choose. 

However, any new explicit exact solutions, even Wronskian solutions, to the KdV 
equation will still be very interesting. Two open questions follow: 

(1) What about the following case? 

JV N 

~<j>i,XX = ^2 /'\.°/.< + ^ 1 - i ~ N > 

3=1 J=l 
N N 

4>i,t = (ij<f>j,xxx + €ij<f>j, i < i < n, 

3=1 3=1 

where all coefficients are real constants. Are there any conditions, rather 
than the ones established in this paper, which will guarantee the Wronskian 
solutions for the KdV equation? 

(2) Are there any similar conditions to guarantee the double Wronskian solu- 
tions for the KdV equation? 

On the other hand, it deserves more investigation whether there exist Wronskian 
solutions for the generalized KdV equations (see, say, [BlIYllifTY] for examples) and 
what differential conditions on Wronskian solutions one can have if there exist 
Wronskian solutions. 
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